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Abstract: The Gaussian curvature of a two-dimensional Riemannian manifold 
is uniquely determined by the choice of the metric. The formulas for comput- 
■^ I ing the curvature in terms of components of the metric, in isothermal coordi- 

Cd ■ nates, involve the Laplacian operator and therefore, the problem of finding a 

Riemannian metric for a given curvature form may be viewed as a potential 
theory problem. This problem has, generally speaking, a multitude of solutions. 
To specify the solution uniquely, we ask that the metric have the mean value 
property for harmonic functions with respect to some given point. This means 
tH- I that we assume that the surface is simply connected and that it has a smooth 

Tij" ■ boundary. In terms of the so-called metric potential, we are looking for a unique 

Cn I smooth solution to a nonhnear fourth order elhptic partial differential equation 

with second order Cauchy data given on the boundary. We find a simple con- 
dition on the curvature form which ensures that there exists a smooth mean 



o 

(•~^ ■ value surface solution. It reads: the curvature form plus half the curvature form 

for the hyperbolic plane (with the same coordinates) should be < 0. The same 

analysis leads to results on the question of whether the canonical divisors in 

^ _ weighted Bergman spaces over the unit disk have extraneous zeros. Numeri- 

H ■ cal work suggests that the above condition on the curvature form is essentially 

^ I sharp. 

• rH . Our problem is in spirit analogous to the classical Minkowski problem, where 

K> I the sphere supphes the chart coordinates via the Gauss map. 

C^ ; AMS Classification: Primary: 53C21, 35J65, 49Q05, 46E22. Secondary: 58E12, 

49Q20, 30C40. 

Keywords: Bordered surface, Riemannian metric, minimal area, mean value 
property, curvature form. 



1 Introduction 

General abstract surfaces. Consider a simply connected C°°-smooth bor- 
dered two-dimensional Riemannian manifold J7. The boundary dfl is then a 



C°°-smooth Jordan curve. We model fl as the unit disk D, supplied with the 
C°°-smooth Riemannian metric ds: 

J7 = (D,(is), ds{z)^ ^ a{z)dx'^ +b{z)dy^ + 2c{z)dxdy, (1.1) 

where we use the convention z = x + iy. The smoothness of the metric means 
that a, b, c are C°°-smooth real-valued functions on the closed disk D, subject 
to the Riemannian metric conditions 

< a{z),b{z), c(z)^ < a{z) b{z). 

We are interested in the problem of reconstructing the metric, if its associated 
curvature form is given. The area form is given by 

(iS(z) = (a(z)6(z) - c(z)2) d^{z), 

where 

dxdy . 

dij(z) ~ , z — X + ly. 

The Gaussian curvature function k is a geometric quantity, which is given by 
Brioschi's formula in terms of the coordinatization as 



(a(0)6(z) - c(z)2)2 1^ 

/ -^d^a{z) - ^dlb{z) + dlyc{z) \dMz) d^c{z) - \dya{z) 
det dyc{z) — ^dxb{z) a{z) c{z) 

\ ia,5(z) c{z) b{z) 

(0 ^dya{z) ^dxb(z) 

^dya{z) a{z) c{z) 

ia,6(z) c{z) b{z) 

The curvature form is given by the expression 

it measures the distribution of the curvature in space. On the curved surface 
ft, there exists a counterpart of the usual Laplacian in the plane, known as the 
Laplace-Beltrami operator, denoted by A. In terms of the given coordinates, it 
can be expressed by 



' dx 



4:y^a{z)b{z)-c{z)^ 



"W ^3.- ,,,.-;« ....a. 



dy 



^a{z)b{z) - c(z)2 ^aiz)b{z) - c(z)2 

c{z) ^ a{z) 



y^a{z)b{z) — c{z)'^ y^a(z)b{z) — c(z) 



We say that a twice diflerentiable function / on the curved surface ft is har- 
monic - or Laplace-Beltrami harmonic, if we want to emphasize that we use 



the Laplacian induced by the metric - provided that ^f{z) = holds through- 
out n. Now, suppose we have two curved surfaces Q, and Q, , which are both 
modelled by the unit disk D: 

rj = (ID,ds), J7' = (D,ds'). 

As the corresponding Laplace-Beltrami operators A and A' are generally dif- 
ferent, we should expect them to give rise to different collections of harmonic 
functions. If they give rise to identical collections of harmonic functions, we say 
that the metrics ds and ds' are isoharmic. This relation between two metrics of 
been isoharmic is an equivalence relation and its equivalence classes are called 
isoharmic classes of metrics. Let 97t denote the collection of all C°° -smooth 
metrics on B, and let dJla run through all the isoharmic classes as a passes 
through a suitably large index set. We then get the disjoint decomposition 

a 

which provides a fibering of Tl. Of particular interest is the fiber TIq (assuming 
that the index set contains the value 0) which contains the Euclidean metric 
as an element. We claim that 97to coincides with the collection of so-called 
isothermal metrics. We recall the standard terminology that ds is isothermal if 

ds(z)2 =: w(z) |rfzp ^Lu{z) {dx^ +dy^), (1.2) 

holds for some C°°-smooth function uj{z) which is positive at each point of D. 
This means in terms of the functions a, 6, c that 

a(z) = b{z) ~ uj{z) and c(z) = 0, 

hold throughout P. The area form is then rfS = cudY^. The Laplace-Beltrami 
operator becomes a slight variation of the ordinary Laplacian 

1 . 



cj(z) 
where A denotes the normaHzed Laplacian 

This means that the Laplace-Beltrami harmonic functions for an isothermal 
metric are just the ordinary harmonic functions. In other words, the isothermal 
metrics are contained in 97lo- To see that all metrics in 9Jlo are isothermal, we 
pick a metric ds in DJIq and note that the functions f{z; zq) = Re [(z — zo)^] and 
g{z; Zq) = Im [(z — zo)^] are ordinary harmonic functions of z, and therefore also 
Laplace-Beltrami harmonic functions. We obtain the equations A.zf{z; zq) — 
and A.zg{z; zq) = and implement the above formula for the Laplace-Beltrami 



operator A. We evaluate the Laplace-Beltrami equation at the point zq; as we 
vary the point zq in D, the assertion that the metric is isothermal follows. 

We would like to better understand the other isoharmic classes OTq. To 
this end, we pick a metric dsi in dJla, for a fixed index a. It is well-known 
in the theory of quasi-conformal mappings that it is possible to find a C°°- 
homeomorphism i^ of D which changes the coordinate chart so that the metric 
dsi becomes isothermal 1 . Under the same coordinate change F, the other 
elements of Tla change as well. However, the collection of Laplace-Beltrami 
harmonic functions, being determined by the geometry of the abstract surface 
in question, remains the same, except for the obvious composition with the 
chart function F. The Laplace-Beltrami harmonic functions for the metric dsi 
are after the coordinate change just the ordinary harmonic functions, and this 
then carries over to all the other metrics in Tla- In other words, the C°°- 
homeomorphism F effects an identification of dJla with 9Jlo . 

The optimization problem. We are given a C°°-smooth real-valued 2-form 
fi on the closed unit disk D, and we are looking for a metric ds on D of the 
type lll.l|l which is smooth up to the boundary and has fi as curvature form 
K = ndTi. There are plenty of such metrics. To reduce their number, we 
decide to minimize the total area of the associated surface J7 — (D, ds) under 
two additional conditions. We require that (1) the area form at the origin is 
essentially the area form of the plane, that is, 

a(0)6(0) - c{Of = 1, 

and we also ask that (2) we only minimize over a fixed isoharmic class Tla of 
metrics. This second requirement means that we fix the collection of harmonic 
functions in the unit disk while performing the minimization. A natural question 
here is whether there exists a minimizing surface, and whether it has a smooth 
boundary. 

In order to solve this problem analytically, it is very helpful to note that 
by using the C°°-homeomorphism F, mentioned in the previous subsection, we 
may choose to work only with the class 9Jlo of isothermal metrics. This means 
that the metrics are of the form 

ds{z)^ =uj{z) |dzp = w(z) (dx^ +dy^), 

for some C°°-smooth function uj{z) which is positive at each point of IJ. The 
formula for the curvature simplifies greatly, 

2 

k(z) = -2Aloguj{z) = --Alogw(z), z G D, 

uj{z) 

and so does the formula for the curvature form: 

K{z) = -2A\oguj{z)d'E{z)^-2A\oguj{z)dJ:{z), z 6 D. (1.3) 



Let us now see how this coordinatization simplifies the formulation of our prob- 
lem. We are given a C°°-smooth real-valued function fi on the closed unit disk, 
which is the density function for the 2-form fi: 

fi{z)^ H{z)dj:{z), zgD. (1.4) 

We write the curvature form K as 

K(z) = K{z) dT.{z), K{z) = -2 A log w(z). 

The equation K = /x then becomes 

-2 A \oguj{z) = ^i{z), z e D. 

As we have already fixed the isoharmic class by setting the Laplace-Beltrami 
harmonic functions equal to the ordinary harmonic functions in the plane, all 
that we really need to fix is the value of uu at the origin, while minimizing the 
total area of the surface Jl = (D, rfs). However, it is equivalent to maximize the 
value w(0) while keeping the total area constant, and we find it convenient to 
fix the latter to equal 1: 

\fl\-s ^ |D|s = / dS(z) = / uj{z) dS(z) = 1. 



The optimization problem and the mean value property. We recall that 
the weight w should solve the following problem: 

maximize tj(0), while 

Aloga;(z) = — - /i(z), z G D, and (OP) 

Lj{z)dJ:{z) = 1. 

Here, /z is a C°°-smooth real- valued function on D, and we ask of to that it too 
should be C°°-smooth (and positive) on D. It is not clear that the optimization 
problem (OP) should have such a nice solution in general. As a matter of fact, 
it is possible to construct counterexamples that are fairly elementary. Never- 
theless, we shall investigate the properties that such an extremal weight u> = loq 
should enjoy. We compare the extremal weight uq with nearby weights tOt, of 
the form 

where h is a C°°-smooth real-valued function on B that is harmonic in the 
interior D, with h{0) — 0, and i is a real parameter. Then 

A\oguJt{z)^ Alogivoiz)^ --^J.{z), zgD, 



and ujt{0) = cjo(O). The extremal property of ujq now forces the inequahty 

wo(z)dS(z)< f uJt{z) dJ:{z) = j e'^^'-^LOo{z)dY.{z) (1.5) 

» Jo Ja 

to hold. By Taylor's formula, 

e*'*(^) = l + f/i(z) + 0(i2), 

for t close to 0. As we plug this into equation ()1.5|l . we arrive at 

LUo{z) d^{z) < / LUo{z) dJ:{z) + t / h{z) ujq{z) dJ:{z) + 0{t^). (1.6) 



By varying t from small positive to small negative values, we realize that the 
only way for lll.(i|l to hold is if 

h{z) a;o(z) dYi(z) = 0. 

If we drop the requirement on h that /i(0) = 0, and consider the function 
h{z) — h{0) instead in the above argument, we obtain 

h{z) loq{z) dT.{z) = /i(0). 

This is what we call the mean value property of the weight ujq. Note that by an 
approximation argument, the above mean value property remains valid when we 
extend the collection of h to all harmonic functions in D that are integrable with 
respect to area measure. We find that we are looking for a (positive) weight loq 
that is C°° -smooth up to the boundary, with 

A\ogujo{z)^--fi(z), zeD, (1.7) 

and the mean value property 

h{z) iOa{z) dY^iz) = /i(0), h e H\B), (1.8) 

where 7i^(D) stands for the Banach space of all complex- valued area integrable 
harmonic functions on D. Strictly speaking, there is more information contained 
in the extremal property of ujq, but clearly, if we find a unique C°° -smooth 
solution Wo to H1.7|l and H1.8|l . then it is definitely our prime candidate for the 
solution to the optimization problem (OP). 

The relationship with Bergman kernel functions. We now discuss how 
to actually find the extremal weight ujq; more precisely, we discuss the problem 
of solving the equations Ijl.Tjl and ljl.811 . First, we note that by elementary 
potential theory, one solution uji to lll.7|l is given by 



logwi(z) = - / log 

Jn 



fi{w)dY.{w), zeB, (1.9) 



and it is well known that logwi is real- valued and C°°-smooth on D, because /i 
has these properties. Any solution to ljl.711 . then, has the form 

logujo{z)^logLUi{z) + H{z), z€D, 

where H is real-valued and harmonic in D; given the smoothness assumptions 
on ujQ, H should be C°°-smooth on B. We find a holomorphic function F on D, 
which is zero-free and C°°-smooth up to the boundary, such that 

log \F {z)\^ = H{z), zeB. 

By restricting lll.8|l to holomorphic functions, we obtain 

/(z)|F(z)|VWdS(z) = /(0), (1.10) 

D 

for all / in ^^(D), the space of area-integrable holomorphic functions on B. Let 
g e ^^(D) be arbitrary, except that ^(0) = 0; then / = g/F is in ^^(D) as well, 
and we find that lll.lfl|l states that 

g{z)F{z)uJi{z)dT,{z)^0. (1.11) 

D 

We interpret this in terms of the Hilbert space ^^(B, wi), consisting of the 
square area-integrable holomorphic functions on D, supplied with the weighted 
norm 

WfU-i I |/(z)pc.i(z)dl](z)" 

This space .4^(D,a;i) is known as a weighted Bergman space. Equation ljl.ll|l 
then states that F is perpendicular to all the functions 

{geA'{D,u;i): 5(0) = O}, 

and this means that F is of the form 

F{z)^CK^,{z,0), 

where C is a complex constant, and Ki^-^{z,w) is the weighted Bergman kernel 
with weight uji . We recall that the weighted Bergman kernel is defined by 

+00 
X^i(z,w) = ^e„(z)e„(w), z,weB, 

where the functions ei(z), 62(2:), 63(2),... run through an orthonormal basis 
for J^{p^u}\) 01 pp. 43-44]. Alternatively, K^^i^-^w) is the unique element of 
^^(D,a;i) which supplies the point evaluation at w S D: 

/(w^) = / Az)K^^{z,w)ujY{z)d^{z), 



for / S -4^(ID', wi) EI- The constant C is easily determined (at least in modulus) 
by applying Ijl.lflll with the choice f — I, which leads to 

F{z) ^ K^,iO,0)-'^^ K^,iz,0), zGD. (1.12) 

Returning back to the extremal weight ujq, we find that it is of the form 

By the elliptic regularity theory for PDEs [Hj, the function i4r„^(-,0) is C°°- 
smooth on ID'. We realize, then, that a necessary condition for the existence of 
a C°° -smooth positive weight loq on D that solves our problem (OP) is that 

K^Jz,0)^0, zgB. (1.14) 

It turns out that it is also sufRcient, and that the extremal solution is then given 
by ^1.1.'^) . For, if w is another weight that solves ljl.711 . and is C°°-smooth and 
positive on D, then it is of the form 

Lu{z)^\F{z)\^LJi{z), ZGD, 

where F is C°°-smooth on fi, analytic in D, and zero-free in D. It is well-known 
that the function 

F{z) = K^,{0,Oy'^^ K^ii^,0), zeD, 

is the unique (up to multiplication by unimodular constants) solution to the 
extremal problem to maximize |i^(0)|, given that F is holomorphic in D and 



\F{z)\^uji{z)dT.{z) ^1, 
which means that cj(0) < cjo(O) for all competitors iv with 

w(z)dS(z) = 1. 



This shows that luq, as given by II1.13|I . is indeed the solution to the extremal 
problem, provided that condition Ijl.l4|l is fulfilled. 

A toy example. We should study a simple example, to develop some intuition. 
We consider the degenerate data 

/i(z) = -9Sx{z), 

where S\{z) stands for the Dirac delta function concentrated at the point A € D, 
and 6* is a real parameter. This does not fulfill our smoothness requirement at 



the point A, so we should think of it as a Hmit case of smooth functions /i„. The 
weight function ui suppUed by equation ljl.9|l is expUcitly given by 



UJi{z) 



A 



1- Az 



z e 



which is a reasonable weight on D provided that —2<9< +00; outside this 
interval, the weight fails to be area-summable near the point A. We need to 
find the weighted Bergman kernel function K^-^ . To this end, we note first that 
if is a Moebius automorphism of D, then we have the following relationship 
between the kernel functions for the weights uj and uj o cj): 

K^o4,{z,w) = (f)' {z) (j)' [w) K^^{4){z) , (t>{w)) , z, w e D. 

Let UJ2 stand for the radial weight 

CJ2(2)=|2|^ zgD; 

then the weighted Bergman space ^^(D, 0^2) has the reproducing kernel function 

1 e 1 

K^^{z,w) = rr^ + TT- r, z,weD. 

(1 — zwy 2 i — zw 

As uji ^ UJ20 <f), where <f) is the involutive Moebius automorphism 

0(z) = -— ^, z e D, 

1 — Az 



we find that 



I e_ l-|Ap 

(1 - zw)2 ^2(1- Az)(l - Aw)(l - zw) '' 



K^^{z,w) ^ — — TT^ + 77 7; T~Tr^ TTTT, ~^' z,we 



Plugging in w = 0, we obtain 

f) \ — I \|2 

K^,{z,0) = l + - U , zeD. 

2 1 — Az 

This expression has a zero in D precisely when (recall that —2<6< +00 is 
assumed) 

In particular, if —2 < 6' < —1, we may choose A close to the unit circle T, 
to make sure that K^-^{z,0) has a zero in D, whereas if — 1 < < +cx3, no 
such zero can be found, no matter how cleverly we try to pick A G B. This 
means that the function uia supplied by relation ljl.1811 is C°°-smooth on ]D)\ {A} 
for —1<0< +CX), which constitutes the solution to our extremal problem. 
However, for —2 < 9 < —1, we see that it is possible to pick A E D so that no 
smooth solution to the extremal problem will exist. 



This calculation suggests a pattern: tor hyperbolic metrics (this means that 
the Gaussian curvature is negative everywhere), we have a smooth solution loq 
to the extremal problem, whereas when the metric becomes elliptic (positive 
Gaussian curvature), we tend to get in trouble. 

Statement of the main results. We are going to compare the curvature form 
with that of the hyperbolic plane, and use this as a criterion for how strongly 
curved our metric is. The Poincare metric is 

, / X 2|dz| 

l-|z|2 

and the associated area form is 

4dS(z) 



d-EMiz) = 



(1-kP) 



2^2 • 



The curvature form for the Poincare metric (which supplies the standard model 
for the hyperbolic plane H) is 

Kh(z) = -- 



(1-|ZP)2' 

we are to compare the curvature form data fji(z) — fi(z) dY.{z) with Kh(z). 

THEOREM 1.1. Let ^ he a C°° -smooth real-valued function on D, and sup- 
pose that the associated 2-form fji{z) = ^(z) d^{z) has 

ti{z) + ^Km{z) < 0, zeD. 

Then the optimization problem (OP) has a unique C°° -smooth positive solution 
070 on B, which is given by relation \l.lfl\) . where uJi is as in \1.0\) . In addition, 

K^,{z,w)j^O, (z,u;) e (DxD) U (DxP). 



This is in line with the intuition we arrived at from our "toy example". Note that 
the assumption of the theorem is considerably weaker than requiring negative 
data /x. As for the necessity of the condition of Theorem II. 11 we have obtained 
the following. 

THEOREM 1.2. Fix a, ao < a < +cx), where uq « 1.04. Consider a C°°- 
smooth real-valued function /i on D, and suppose that the associated 2-form 

/x(z) — ^,(z)dYi(z) has 

/x(z) + -Kh(z) <o, zero. 

Then there exists a choice of /i such that the optimization problem (OP) fails to 
have a C°° -smooth positive solution loq on B. 



10 



We conjecture that Theorem II .21 will remain true with ao = 1, making the 
statement of Theorem II . II essentially sharp. 

Work related to the problems considered here can be found in the papers 

i, 0, 0, 0, [HI, [EI, m, [m, ini. 

The metric potential. The Green function for the Laplacian A is 

2 



G{z,w) = log 



z — w 



z,w £ 



z ^ w. 



1 — zw 
The metric potential associated with the isothermal metric ljl.211 is the function 



^{z)^ G{z,w)Lu{w)dJ:{w), zeP; 

see, for instance, ^^I- It solves the boundary value problem (T is the unit circle) 



A#(z) = 1, z e . 

*(z) = 0, z G T. 



(1.15) 



It is set in boldface because it expresses a quantity that is independent of the 
choice of coordinates. We wish to describe the potential equation ljl.711 and 
the mean value property lll.8|l . which the smooth solution to the optimization 
problem (OP) should satisfy, in terms of the metric potential. If we let $o stand 
for the metric potential associated with a weight ojq with H1.7|l and H1.8II . we 
find that $o solves 

A\ogA^a{z) = --fi{z), 



*o(2) = 
d 



2' 

zeT, 



z e 



(1.16) 



*o(^) = 2, z G T, 



dn{z) 

where d/dn{z) denotes the normal derivative, taken in the exterior direction. 
We are of course only looking for subharmonic solutions $0; which means that 
the expression log A$o is more or less well-defined (at least if we have some ad- 
ditional smoothness, and the subharmonicity is "strong"). The normal derivative 
condition in lll.lfi|l cleverly encodes the mean value property H1.8|l . as is seen 
easily from an application of Green's formula. We realize that Theorem 11.11 
can he interpreted as an assertion claiming the existence of a unique smooth 
solution to the non-linear elliptic boundary value problem ^.16)) under appro- 
priate conditions on ^. By moving the point at the origin around by applying 
a Moebius transformation that fixes the unit disk, we find that the method can 
also treat the case when the data for the normal derivative is replaced by 



d 



*o(^) 



1-|A|= 



dn{z) "' ' |A-z|2' 

where A is any point of D. It would be interesting to have an analysis of the 
equation Ijl.l6|l . where the boundary data are considered from a wider class of 
functions. 
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2 Preliminaries 

The standard weighted Bergman spaces. For 1 < a < +00, let ^q(ID') 
A'^{^,uJa) denote the weighted Bergman space for the weight 

w„(z) = (a + l)(l-|zn", zeD. 

The norm in A^ (D) is written 

1/2 






Take a point A £ D \ {0}, and let (p\{z) — ip\,a{z) be the so-called extremal 
function for the problem 



sup 



{Re/(0): /(A) = 0, ||/|U < l}, 



which is unique, by elementary Hilbert space theory. If we let K\ = K'^ denote 
the reproducing kernel for the closed subspace 

{feAlm-. /(A) = o}, 

then the extremal function can be written 

ifxiz) - Kx{0, 0)-i/2 ^^(^^ 0), z e D. 

When the function ip\ serves as a good divisor of the zero at A in the space 
^^(D), it is called the canonical divisor oiX as it is called in [S], or the contractive 
zero divisor, provided that division by it defines a norm contractive operation 
as it is called in |2j. It is known ^| pp. 58] (see Section 4 for details) that the 
reproducing kernel K — Ka for the space A^ (D) has the form 

(1 — zw)"+^ 
and that the kernel K\ is derived from K via the identity [T?i| 

K{z,\)K{\w) 



K\{z, w) = K{z, w) 



1 (i-iAn"+2 



(l-zw)"+2 (l-zA)"+2(l-Au;)"+2' 
(see Section 4 for the details) . It follows that the extremal function is given by 

^.(z) . (1 - (1 - \x\-rr"'^ {1 - ^^37^} ' - ^ ^- (2-1) 

The following property of the extremal function Lp\ — tpx^a is fundamental. 
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LEMMA 2.1. Fix a in the interval — 1 < a < +00. Then, for each bounded 
harmonic function h on D, we have that 



h(z) i^A WP (« + 1) (1 - kP)" rfsw = MO). 



Proof. First, suppose h is an analytic polynomial. Then, if we recall the defini- 
tion of ip\ in terms of Kx , and use the reproducing property of the kernel K\ , 
we have that 



h{z)\^x{z)Wa + l){l~\z\^Yd^{z) 

M^)|^^a(z,0) (a + 1)(1 - |zn"dS(z) = MO). (2.2) 

Taking complex conjugates in Ij2.2|i . we realize that the desired equality holds 
for all harmonic polynomials (defined to be sums of analytic and antianalytic 
polynomials). A simple approximation argument finishes the proof. D 



A weighted biharmonic Green function. The biharmonic Green function 
is the function F on D x D that solves the boundary value problem 



A^r(z,u;) = 5w{z), z G 



r(z, w) — 0, 

d 



^ dn{z) 



zeT, 
r(z,w) = o, zeT, 



(2.3) 



where to G D, and 5w stands for the unit point mass at the point w. We will 
think of locally summable functions / on some domain Q, of the complex plane 
as distributions on 17 via the linear duality 



(/» 



/(z)0(z)dl](z). 



for compactly supported test functions cf) on O. Given this normalization, the 
biharmonic Green function is given explicitly by the formula 



r(z, w) — \z — w\ log 



z — w 



1 — zw 



+ (i-izn(i-i^n, 



z,w G 



We shall also need the Green function for the weighted biharmonic operator 

A(l — |zp)^^A, denoted Fi, which, by definition, solves, for fixed w G B, 



f A{1 -\z\Y^ AT, {z,w)^5^{z), 

Ti{z, w) — 0, z G T, 

d 

ri(z, w) =0, z G T, 



z G . 



(2.4) 



. dn{z) 
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Although the differential operator is singular at the boundary, the above bound- 
ary value problem has a unique solution (we may use Green's theorem to inter- 
pret the boundary data in terms of integral conditions for AzTi{z,w), which 
are uniquely solvable). 

The function Fi was calculated explicitly in [2]: 

LEMMA 2.2. We have that 



ri(z,u;) = <j |z - wp - -|z^ - w^l [►log 



1 



z — w 



1 — zw 



+ - (1 - \z\')il - \w\')i 7 - \z\' - \w\' - \zw\' - 4Re ziu 



2(l-|zn(l-H^)| '""'I' 



|1-ZW|2 



for z, w e D. Moreover, it follows that, for z, w e D, 

1 (l-|z|2)3(l_|u.|2)3 



|l-ziD|2 



< 



■ <Ti{z,w) 

1 ii-\z\^rii-\w\'r 



8 



|l-zw|4 



||l-zwp+4- |z + ?«p|. 



Proof. Let Fi^stand for the function defined by the above expression; we want 
to show that Fi = Fi. To this end, we show that it solves the boundary value 
problem which determines Fi. We first note that 

A^ri(z,'u;) = {I - \z\^)[G{z,w) + Hi{z,w)], z,w £D, 

where Hi{z, w) is the harmonic function of z which is given by 

1 .- , ,^^ 1 — |zWp ,^ , ,9s „ r zw 



H,iz,w) = {l-\wn\-i3~\w\') ' _ +(l-|u;nRe 

1^2 |1 — zw|^ 1(1 — zw)^ 

for z,w gD. Thus, Fi satisfies 

A(l - |z|2)-iAfi(z, w) = S^{z), z e D. 

We shah now estabhsh the specified inequality for Fi, which entails that 

fi(z,^) = 0((l-|z|)3), |z|^l-. 



It fohows that for a fixed w G B, the function Fi satisfies the boundary condi- 
tions 

f Fi(z,w) =0, zeT, 

I 9„(2)Fi(z,w) =0, z e T, 
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where dn{z) is the outer normal derivative. In view of this, we conclude that 
Fi = Fi. It remains to establish the claimed bounds for Fi, from above and 
from below. We use the following estimate of the logarithm: 



r 1 3 r^ 



0<r < 1, 



with the choice 



1 — zw 
and base our calculations on the identity 



1 



z — w 



1 — zw 



(i-izn(i-H^) ^ 

|1 — zw\^ 



The steps are rather lengthy but quite elementary, and are therefore left to the 
reader as an exercise. D 

The next result is a consequence of the positivity of the Green function Fi. 

LEMMA 2.3. Let u be a C^ -smooth suhharmonic function on B such that for 
some real (3, < (3 <2, 



Hz)\ = o 



1 



(i-kl 



1" 



Also, let ip\ = ipx^i denote the extremal function for the point A G D m the space 
Al(B). Then we have the following inequality: 

u{z)il - \z\^)d^{z) < f \ipx{z)\^u{z){l ~ \zf)dl^{z). 

Proof. We consider the potential function 

<Pxiz) = f G{z,w) (Iv^aMP - 1)(1 - kHriSH, z e B, 



which solves the problem 

r A$,(z) = (|^,(z)p-l)(l-|z|2), zeD, 
j $aW = 0, zeT. 

Note that by using Green's formula as in ,7"., we see that the property that (p\ 
has according to Lemma lTTl mav be rephrased in terms of the potential function 
$a: 

-4-^$a(2) = 0, zeT. 
an[zj 
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It follows that $A solves the over-determined boundary value problem 

A<&a(z) = (|^a(z)|2-1)(1-|z|2), ;^ep, 

$a(z) = 0, z e T, 

We may remove this over-determination by increasing the degree of the elliptic 
operator: 



A^ —A^xiz)^\^'^iz)\^ zeB, 



1 
d 



. dn{z) 



$a(2) = 0, zeT. 



This problem has a unique solution, which may be expressed in terms of the 
Green function Fi, 

<Px{z) = f ri(z,w) |(^^H|2dEH > 0, z e B. 

Jn 

We first apply this to the case when u is C^-smooth on D, and see that in view 
of the assumption that u is subharmonic, we find, by an application of Green's 
theorem, that 

{\^x{z)\^ - 1) u{z) (1 - Izp) dS(z) - / $aW Au{z) dnz) > 0. 

If u is not smooth up to the boundary, we perform the above for the dilated 
function Ur{z) = u{rz), with < r < 1. By the growth assumption on u and 
the smoothness of the function ip up to the boundary, we may let r ^ 1~ and 
apply Lebesgue's dominated convergence theorem, to conclude that 

{Wx{z)?-l)u{z){l-\z\^)dnz)>Q. 

The proof is complete. D 

Now, let w be a strictly positive C°°-smooth weight on D, which is such that 

the function 

uj{z) 

^ ^ i°s 'a — ri2 

1- |z|2 

is subharmonic on D. As before, let K^^(z,w) denote the reproducing kernel 
for the weighted Bergman space J?(Ii,uj). The following result is basic for our 
further considerations. 

THEOREM 2.4. The function K^ extends to he C°° -smooth on the set (© x 

D)U(P xP). 
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This is a classical theorem of elliptic regularity type, obtained in 1955 by L. 
Nirenberg ^j. It is independent of the above subharmonicity requirement. 

Next, we consider the function A^^ : D -^ K, as given by 

|g^(z,0)P u:{z) 

it is positive and subharmonic in D, due to the assumptions on w, and it has 
the growth behavior 

A basic property of A^^ is the following. 

LEMMA 2.5. For each hounded harmonic function h on D, we have that 



h{z) A^(z) 2(1 - |zp) dS(z) = h{0). 

Proof. First, let us assume that h is an analytic polynomial. Then, using the 
reproducing property of the kernel function K^^ , we have that 

/i(z)A^(z)2(l-|z|2)dS(z) 

/^(^)l^ll^^^(^, O)^(z) dll{z) = h{0). (2.5) 

Taking complex conjugates in (j2.5|l . we see that the desired equality holds for 
all harmonic polynomials. An approximation argument finishes the proof. D 

This means that the function A^^ has a lot in common with the function \(p\\'^ 
for the parameter a = 1, which suggests it may have an expansive multiplier 
property that is similar to the one obtained in Lemma [2. 31 

LEMMA 2.6. Suppose that u is a C^ -smooth subharmonic function in D, 
which has the growth bound 

I"(^)I^q( (i_I^I)/3 J «« 1^1 ^1"- 

for some real /3, < /3 < 1 . We then have 

u{z) 2 (1 - Izp) dS(z) < / A^{z)u{z) 2 (1 - |zp) dE{z). 
) Jb 

Proof. We introduce the potential function 

$^(z) = / G(z, w) [A^{w) - 1] 2 (1 - \w\^) d^{w) 

Jo 



17 



which solves the boundary value problem 

f A$^(z) = 2(l-|z|2)[A^(z)-l], zeB, 

\ $^(z)-0, zeT. 

It follows from Green's formula and Lemma ITKl as in 0, that the potential 
function has 

---$„(z) = 0, zeT. 

on(z) 

Then ^^ solves the over-determinated boundary value problem 

A$^(z) = 2(l-|z|2)[A^(z)-l], ze© 

$^(z) - 0, z e T, 

-4-^$^(z) - 0, z G T. 

Increasing the degree of the elliptic operator, we find that <i>cj also solves the 
problem 

A(l-|z|2)-iA$^(z)==2AA^(z) >0, zeD 

$^(z) ==0, z e T, 



9n(z) 



$^(z) = 0, z e T, 



which has a unique solution. Then ^ui may be expressed in terms of the Green 
function Fi, 

$^(z) = 2 / ri(z, w)AA^(u;) dS(u;) > 0, z e B. 

Let us consider now that u is a subharmonic function which is C^-smooth on 
B. Then, by applying Green's theorem, it follows that 

(A^(z) - 1) u(z) 2 (1 - Izp) dS(z) = / $^(z)Aw(z) dE(z) > 0. 



In the case when u is not smooth up to the boundary, we consider the dilated 
function Ur{z) = u{rz), with < r < 1, for which the above inequality holds. 
It follows from the growth bounds of u and A, and Lebesguet's dominated con- 
vergence theorem, that we may let r ^ 1~ to conclude that 



(A^(z) - 1) u{z) 2 (1 - \z\') d^{z) > 0. 
The proof is complete. D 
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3 The proof of Theorem 11.11 

We realized back in the introduction that, in order to obtain Theorem ll.il all 
we need to do is show that 

K^{z,w)^0, (z,w) e (D xD) U (D xB), 

provided that w is a C°°-smooth and positive weight function on D, with the 
property that 

uj{z) 

^ ^ i°g 'i — ri2 

1- |z|2 

is subharmonic on D. After all, Theorem 12.41 guarantees that the weighted 
Bergman kernel is C°° -smooth up to the boundary. It is easy to verify that the 
above subharmonicity requirement is the same as the condition on the curvature 
form in the statement of Theorem ll.il 

The proof splits naturally into two parts. 

PROPOSITION 3.1. Under the above conditions on lu, 

K^{z,w):^0, {z,w)eBxB. 

Proof. We observe first that for any Moebius map preserving the disk D, we 
have 

K^{(l){z),(j>{w)) ^ K^^iz,w), z,we]D), 

where 

LU^{z)^\(t)'{z)\'^LUO(l){z), ZeD. 

We claim that oj^ is a weight of the same type as w. In fact, the function 

z h-^ log ' 



is subharmonic on D if, and only if, the function 

uj{z) 



log 



,1-|0-1(Z)|2, 

is subharmonic on D as well. Then, if we consider a Moebius map 

r'W = 7f^, 2,Cee, |7| = 1, 

1- z( 



we find that 






,i-i0-i(2)iv "vi-kiv "Vi-ici 

Thus, the function 

z H^ log ' 



\-\z\ 
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is subharmonic on D if, and only if, 



log 



is subharmonic on D as well. It follows that it is enough to specialize to w = 0: 

if^(z,0)^0, zgE). 

We note that, by the reproducing property of the kernel function, ^^^(O, 0) — 
\\K^{z, 0)1]"^ > 0. We introduce the extremal function L, given by 

L{z) = {K^{0,0))-^K^{z,0), zeD, 

which solves the problem 

sup{Re/(0): ||/|U < 1}. 

By Theorem ITil the function L is C°°-smooth on D. 

We argue by contradiction. So, we assume that there exists a A G B such 
that K^{X,0) — 0; then L{\) = 0. Consider then the function 

L{z)^L{z)/^xiz), zeD 

where ipx is the canonical divisor of {A} in the space ^^(ll]'). We should point 
out that, due to the smoothness of ip\ and the fact that it doesn't have any 
extraneous zeros on D, the function L is also C°°-smooth on D. 
Let u be the function given by 



2(l-|zP)l 



<^) = :T7P^Tl^\Liz)\', ^ ^' 



It follows from the hypothesis on the weight that log m is a subharmonic function 
in D, so that in particular, u is subharmonic as well and it has the growth bound 

It follows from Lemma, l273l that 

||L||^= I \L{z)\^u{z)dnz) 

u{z) 2 (1 - Izp) d^{z) < f |(^A WP u{z) 2 (1 - Izp) dE(z) 

\Liz)\^u;{z)d^iz)^\\L\\l^l. 

On the other hand, it follows from equation Ij2.1|l that ip\{0) < 1 and so 
1/(0) > L{0), which violates the extremal property of L{z). This is the de- 
sired contradiction. Hence, the function L does not have zeroes in D. D 
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PROPOSITION 3.2. Under the above conditions on u, 

K^{z,w)^0, (z,w) e (T X D) U (D xT). 

Proof. Note that, by the same argument used in the proof of the Pronosition l.S.ll 
and the fact that 

K^{z,w) ^ K^{w,z), z,weD, 

it is enough to prove that 

K^iz,0)^0, zeT. 

We argue by contradiction. So, we shall assume that there exists A G T such 
that Kuj{X,0) — 0. Then, due to the smoothness of z i-^ Ku;{z,0) on fj), as 
provided by Theorem 12.41 we find that 

\K^{z,0)\=O{\z-X\), as D3z->A. 

Let < r < 1 and define the function fr G ^^(ID)), given by 

^ ' ^Ki{rX,rX) (l-rAz)3' 

where Ki is the reproducing kernel for the space AKB). It follows that 

||M|i = l, 0<r<l, 

where || • ||i is the norm in Ai{lS>), as defined back in Section 2. Furthermore, 
l/rp is bounded on B for each < r < 1. We now consider the function 

It follows that 

Ru^iz) = 0{\z-Xf) as B3z^X. 

Then there exists a positive constant M such that 

R^{z) <M\z- Xf, zGD. 

It follows from Lemma [2. 61 and the inequality 

r-jz- A| < |l -rAz|, z G D, < r < 1, 

that 

1= f \friz)\' 2(1- \zf)dJ:{z)< f AM)\fr{z)\'2il-\z\')dJ:iz) 



= /"i?^(z)|/,(z)pdS(z)<M /|z-An/,(z)|2dS(z) 

Jvi Jo 



(1 - r^)"" f 1 , , , 1 - r2 

|l-rAz|4 



< M '--^ I ^-=-^dS(z) = Af „ 



which is a contradiction for r sufficiently close to 1. D 
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4 The proof of Theorem 11.21 

In this section, we construct, for a fixed a > ao — 1.04, an explicit example of 
a function ^ € C°°(D), with the associated 2-form /i — /irfE, such that 

/x(z) + |Kh(z)<0, zeB, (4.1) 

for which the optimization problem (OP) fails to have a smooth-positive solution 
Wo on D, as the weighted Bergman kernel Ki^^{-,{)) has an extraneous zero in 
ID). Here, ui is associated with /i as in the introduction: 



\oguji{z) = - / log 



1 — zw 



fiiw)dT.{w), zeB. (4.2) 



The choice of /x. We first consider the extremal case for the inequality IJ4.1II . 



a 2 a (iE(z) 

2^«('^=(l-|z|2)2 



Mh(2) = -^7 Kh(2) = p^, z G 



In this case we can compute explicitly the weight wq , which solves the optimiza- 
tion problem (OP), 

t^o(^) = (1 - kP)", zgB. 

Following the same line of thought as in our toy example from back in the 
introduction, we then consider the data function 

where A = {ak}k is a finite collection of points in B, {pk}k is a convenient 
sequence of positive constants, and Sa^. stands for the Dirac delta function, 
concentrated at the point a/c e B. The 2-form fj, — fidH meets the inequality 
H4.1II , but fi is very rough at the points of A and hence it does not satisfy the 
C°°-smoothness requirement. However, it is easy to approximate a point-mass 
by a sequence of positive C°°-smooth functions. Moreover, we may replace the 
function 

''^^'~ (l-|zP)2 

by a slight dilation. 



Vr{z) 



2r^ a 



(l-r2|0|2)2 

with r, < r < 1, close to 1. This means that if 



K^,{z,w)^Q, z,u;eB, (4.4) 
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holds for smooth indata ^, with fi and uji connected via IJ4.2II . then it also 
holds for rough indata of the above type, modulo some slight modifications. To 
explain these modifications, we note that 



log 



1 — zw 



i>j.{w) dY,{w) =alog 



1-H 

l-r^\z[ 



z G. 



Let the weight wi.r be defined by the formula 114. 211 . where iO\ is replaced by 
wi^r-, and [i is replaced by 



fj.r{z) 



2r'^a 



(l-r2|z|2)2 



Xl^fc^'»'=(^)' 



z G 



We then calculate that 



u;i,(z) = (l-rr"(l-r^N^)"n 



z - ak 



1 - QkZ 



Pk 



z G. 



As reproducing kernel functions have homogeneity index —1 in general, that is, 

Ktu,iz,w) = -K^{z,w), 

holds for arbitrary positive constant t, we see that 

K^,Jz,w)^{l-r^rK^,Jz,w), 
where uJ2,r is the weight 



u;2A^) = {l-r'\z\Y'[[ 



z - ak 



1 - QkZ 



Pk 



z G . 



If 114. 4|l holds for uji = uji^r, then it also holds for wi — W2,r, and vice versa. As 
r ^ 1~ , the weight LU2,r tends to 



u;2{z) = {I - \z\ril 



z - ak 



1 - akz 



z G 



(4.5) 



and if the reproducing kernel functions for the weights uj2,r are all zero-free 
in ro^, then so is the reproducing kernel for the weight uj2, by a limit process 
argument. We shall prove that with appropriate choices of the configuration of 
the points A — {ak}k as well as of the positive parameters pk, the reproducing 
kernel function for uj2 will have zeros in D^. This then shows that also with 
smooth data, we must have zeros in the associated reproducing kernel function. 
We should note that uj2 satisfies 

A log UJ2 {z) = --fi{z), zgB, 



23 



where /i is given by H.l^ . This means that uj2 is the weight we where looking 
for. 

For computational reasons, we shall only consider pk such that pfe/2 is a 
positive integer. 

Relations between kernel functions. A closed subspace / oi A^{I}) is called 
invariant if it is invariant under the multiplication by the identity function, more 
precisely, ii zf £ I whenever f G I. For a sequence A = {ai, 02, . . .} of points 
on D, the subspace Ia consisting of all functions in ^^(D) whose zero sets 
contain A, counting multiplicities, is an invariant subspace. Such a subspace Ia 
is called zero-based invariant subspace. The reproducing kernel function for the 
invariant subspace Ia is as usual defined by the formula 

+00 
K%{z,w) = ^en{z)en{w), z,weB, 

n=0 

where the functions ei(z), 62(2;), 63(2), . . . form an orthonormal basis for Ia- It 
has the reproducing property 

/(z)= f K^{z,w)f{w)ia + l){l-\w\Ydnw), zeB, 

for f (z Ia- For a finite subset A — {ak}k of D, the associated (finite) Blaschke 
product is the function 



Ba{z) = 1[^^, ze 

^■^ 1 - ahz 



k 

We consider the following two weights: 

u^iz) - (a + 1)(1 - |z|2)", uj^^a{z) - (a + 1)(1 - \z\^)" \Ba{z)\\ 

Note that ijJa,A equals the weight W2 as defined by H4.5II . with all the parameters 
Pk set equal to 2. 

The following proposition is well known. 

PROPOSITION 4.1. We have the following identity of kernels: 

K%{z, w) = Ba{z) Ba{w) K^^,Az, w), z,weB. 

In view of the above proposition, we need to look for extraneous zeros in the 
reproducing kernel function for Ia in order to get zeros of the kernel function 
for the weight UJ2 = u}a,A- 

The kernel function for a zero-based invariant subspace. When the 
sequence A consist of a finite number of distinct points, the kernel function K')^ 
for I A may be obtained by means of the well-known iterative formula (see [9j 
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The first step of this iteration is to apply the formula 114. 7|l to the case A = 9 
and A £ D, to get 

1 n_U|2)a+2 

In the case where the zero set A contains repeated points, i.e. zeros with high 
multiplicity, the kernel function K^ for Ia can be computed via an iterative 
formula, similar to H4.7II . but involving the derivatives of the kernel. Namely, if 
we assume A has no multiple points, then 

^0°(^'^)^ (1 - zzIJ)("+^) ' ^'^e®' 

azUwJ^A\z, w)\^^^^ ^^^ 
where 

When A consist of n copies of the point a G D, the reproducing kernel function 
K^ for the subspace Ia is 



K%{z,w)= ^ 



(1 - zu;)"+2 



n— i . _ _ . 

E / z — a Y / w — a Y 



where 

_ r(a + 2 + j) 

Our computational work suggest that whenever a is bigger than 3, it is pos- 
sible to find extraneous zeros; for a < 3, however, it seems that such extraneous 
zeros do not occur. 

Computational implementation. The recursive formula for 114. 7|l is some- 
what inconvenient in computational applications. Instead, we use the fact that 
the kernel function with w = may be expressed as 

^3(^^0) = 1-^(13^' (4.8) 
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where the {cj}j are certain complex coefficients, which depend on the sequenced 
and on the parameter a. To determine them, we need to solve the linear system 
K2{aj,0) = 0, for every aj G A. This we rewrite in terms of a square matrix 
M, whose entries are given by 



^''^ = n -^-77.^^+2 ' °^' °J ^ ^' 



1 

the equation that determines the coefficients is the linear system 

Mc = l, (4.9) 

where c = {cj}j is the coefficient vector in column form, and 1 is the column 
vector with the number 1 in all its coordinates. 

Although the numerical computation of the kernel function K^ for several 
zeros via the linear system 114. 9|l is technically possible, we are to treat a rather 
badly scaled matrix. Numerically, the matrix M is nearly singular, and this 
leads to possibly large numerical error while solving the system Ij4.9|l . In ad- 
dition, the kernel function expressed as expressed by (j4.8|l is quite sensitive to 
even small perturbations in the appearing coefficients already when the number 
of zeros in A is rather modest. Taken together, this forces us to work with 
higher precision than what is standard. We use 40 decimal digits of precision, 
instead of the standard 16. Most of the numerical tests were done in MATLAB 
6.1. To work with forty decimal digits, we used the command VPA (variable pre- 
cision arithmetic). We estimate the computational error in the approximated 
kernel function K^i-, 0) by comparing the computed values at the points of the 
given zero set A — {ak}k, where the function vanishes. If these values are suf- 
ficiently small, we may be certain that the negativity of K'^{1,0) is a genuine 
phenomenon, provided that the negative value is substantially bigger than the 
numerically obtained values at the given collection of zeros. 

The appearance of an extraneous zero. We first observe that if the ele- 
ments of the zero set A are distributed symmetrically respect to the real axis, 
then the kernel function K^ix^O) is real-valued for real x. Since, trivially, 
< K2{0,0), and since the function K^iz^O) is continuous in the closed unit 
disk D, it follows from the Mean Value Theorem of Calculus that K'^{xo,0) = 
holds for some a;o, < a;o < 1, provided that i^3(l, 0) < 0. 

After testing several patterns for the distribution of the zero set A (including, 
for instance, a multiple zero at a single point), we focused our computations on 
a configuration which yields extraneous zeros for a all the way down to 1.04. 
The same pattern seemed to emerge also when the computer was allowed to pick 
the configuration of the given zeros according to a so-called genetic algorithm. 
The configuration depends on the number n of points of A, as well as on two 
parameters 9 and d, with < 9 < ^n and 1 < d < +oo. By the construction, n 
is an even number. The points of A are given by 

ak = exp {3(i - 9) d'=-"/2}, fc = 1, 2, . . . , n/2, 

ak+^ = CLk, k = 1,2,. ..,n/2. 
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Figure 1: Zero set A for n = 80 generated by the parameters 9 — 0.15, d — 1.30, 
for a — 1.50. The indicated straight lines are tangents to the curves on which 
the zeros are located. 



Figure 2: (top) Level curves of the kernel function iir^(-,0) for the zero set 
A generated by the parameter values: a = 3, n = 6, — 0.51, d = 10. 
(bottom) Values of the kernel function K'^{1, 0) for the zero set A generated by 
the parameters: n = 1500, 6 = 0.033, d = 1.045 and 1.035 <a< 1.050 



In Figure Q we plot the configuration of the given zero set in the complex 
plane for n = 80 zeros, with the parameter values a = 1.5, 9 ~ 0.15, and 
d= 1.3. 

In Figure [2(top), we plot the level curves around the extraneous zero of the 
modulus of the kernel function for the zero set A generated by the parameter 
values a — 3, n = 6, 9 = 0.51, and d = 10. In Figure |2l(bottom) , we plot the 
kernel function K'^{1, 0), while varying a in the interval [1.035, 1.050]. 

Inspiration for our numerical work was derived from ^2]) where Jakobsson 
obtained extraneous zeros for a « 1.40. 
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In the table below (Table 1), we supply values of n, 0, and d, for which our 
computations indicated that the kernel function possesses extraneous zeros for 
a prescribed value of the parameter a. 



a 


n 


9 


d 


3 


6 


0.51 


10 


2.5 


8 


0.48 


8 


2 


14 


0.351 


3 


1.6 


26 


0.265 


2.1 


1.25 


78 


0.176 


1.52 


1.118 


230 


0.104 


1.22 


1.1072 


272 


0.092 


1.183 


1.097 


340 


0.07725 


1.141 


1.065 


550 


0.07 


1.13 


1.053 


770 


0.0556 


1.09 


1.046 


944 


0.0497 


1.078 


1.043 


1090 


0.0445 


1.067 


1.04 


1500 


0.033 


1.045 



Table 1. 



Hinted analytical solution of the problem. The configuration of zeros 
in the above numerically-based counterexamples which works at least down to 
the parameter value a w 1.04 suggests that analytically, we should "smear out" 
the zeros along the two lines; we may decide to work almost infinitesimally 
close to to the point 1, and by blowing up, we may assume that the domain 
is the upper half plane, where the origin plays the role of the point 1. We 
place hyperbolically equi-distributed smeared-out "zeros" along two half-lines 
emanating from the origin, symmetrically located with respect to the imaginary 
axis. In geometric terms, this means that we construct a new weight - the old 
one being (Imz)" - which is the same as the old weight in the region between 
the two half-Hnes, but is the old weight times the exponential of a constant times 
the angle to the nearest half-line in the remaining two regions near the real line. 
We suspect that the reproducing kernel for this new weight (with one argument 
fixed equal to i) has a zero somewhere along the imaginary axis for each fixed 
a, 1 < a < +0O, provided the angle of the two half-lines to the real fine are 
chosen appropriately, and the constant that regulates the density of the "zeros" 
is appropriate as well. However, the actual implementation of this scheme is 
not easy, because it is generally a hard problem to calculate reproducing kernel 
functions for weights that do not exhibit strong symmetry properties. We would 
Hke to be able to return to this problem. 
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